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Abstract. Cosmological implications of clusters of galaxies are dis- 
cussed with particular attention to their importance in probing the cos- 
mological parameters. More specifically we compute the number counts 
of clusters of galaxies, LogiV-LogS' relation, in X-ray and submm bands 
on the basis of the Press-Schechter theory. We pay particular attention 
' to a set of theoretical models which well reproduce the ROSAT 0.5-2 

keV band LogiV-LogS' , and explore possibilities to break the degeneracy 
among the viable cosmological models. 
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1. Introduction 



O ' There are several reasons why clusters of galaxies are regarded as useful probes 

. of cosmology including (i) since dynamical time-scale of clusters is comparable 

^ | to the age of the universe, they should retain the cosmological initial condition 

fairly faithfully, (ii) clusters can be observed in various bands including optical, 
X-ray, radio, mm and submm bands, and in fact recent and future big projects 
^ ■ (e.g., SDSS, AXAF, PLANCK) aim to make extensive surveys and detailed 

imaging/spectroscopic observations of clusters, (iii) to the first order, clusters 
are well approximated as a system of dark matter, gas and galaxies, and thus 
theoretically well-defined and relatively well-understood, at least compared with 
galaxies themselves, and (iv) on average one can observe a higher-z universe with 
clusters than with galaxies. In particular X-ray observations are well-suited for 
the study of clusters since the X-ray emissivity is proportional to n\ and thus 
less sensitive to the projection contamination which has been known to be a 
serious problem in their identifications with the optical data. 

In fact, various statistics related to the abundances of clusters has been 
extensively studied to constrain theories of structure formation, including mass 
function (Bahcall & Cen 1993; Ueda, Itoh, & Suto 1993), velocity function (Shi- 
masaku 1993; Ueda, Shimasaku, Suginohara, & Suto 1994), X-ray Temperature 
function (hereafter XTF, Henry & Arnaud 1991; White, Efstathiou, & Frenk 
1993; Kitayama & Suto 1996 ; Viana & Liddle 1996; Eke, Cole, & Frenk 1996; 
Pen 1996). Previous authors have focused on cosmological implications of clus- 
ter XTF mainly because theoretical predictions are relatively easier although 
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Figure 1. Left: Theoretical predictions for LogiV-LogS' of X-ray 
clusters in CDM models with different cosmological parameters; (a) 
as = 1-04 models with different Qo, Ao and h, (b) Oo = 1 and 0.45 
models with different a$. Denoted by (COBE) are the models normal- 
ized according to the COBE 4 year data (Bunn & White 1997). Data 
points with error bars at S £ 1CP 12 erg cm~ 2 s _1 are from the ROSAT 
Deep Cluster Survey (RDCS, Rosati et al. 1995; Rosati & Delia Ceca 
1997), and the error box at S ^ 2 x 10~ 12 represents a power-law fitted 
region from the R OSA T Brightest Cluster Sample (BCS, Ebeling et al. 
1997). For the BCS data at S = 2 x 1(T 12 , 1 x KT 11 and 6 x KT n erg 
cm -2 s , we also plot the corresponding Poisson errors. Right: Limits 
on and as in CDM models (n = 1, h = 0.7) with (a) Ao = 1 — ^o, 
and (b) Ao = 0. Constraints from cluster LogA^-LogS (solid) and XTF 
(dotted) are plotted as contours at 1ct(68%), 2ct(95%) and 3<r(99.7%) 
confidence levels. Dashed lines indicate the COBE 4 year results from 
Bunn & White (1997). 

the observational data are statistically limited. In addition, the conversion to 
the number density at high z becomes very sensitive to the adopted cosmological 
parameters. On the other hand, LogiV-LogS' which we discuss in details below 
is observationally more robust (Ebeling et al. 1997; Rosati & Delia Ceca 1997) 
while its theoretical prediction is more model-dependent (Oukbir, Bartlett, & 
Blanchard 1996; Kitayama & Suto 1997; Kitayama, Sasaki & Suto 1998). In 
this respect, both statistics are complementary. 



2. Log N — Log S of X-ray clusters 



We compute the number of clusters observed per unit solid angle with X-ray 
flux greater than S by 



N(> S) 



dz d\ (z) c 



dt 

dz 



/ dS {l + z) 3 n M (M,z) 
Js 



dM dT gas <iL band 



dT gas dL 



dS 



where c is the speed of light, t is the cosmic time, d A is the angular diameter dis- 
tance, T gas and L band are respectively the gas temperature and the band-limited 
absolute luminosity of clusters, and hm(M, z)dM is the comoving number den- 
sity of virialized clusters of mass M ~ M + dM at redshift z. 

Given the observed flux S in an X-ray energy band [E a ,Ef,], the source 
luminosity L band at z in the corresponding band [E a (l + z),Eh(l + z)\ is written 
as 

L band [E a (l + z),E b (l + z)}=47Td L 2 (z)S[E a ,E b ], (2) 

where d L = (1 + z) 2 d A is the luminosity distance. We adopt the observed 
L bol — T gas relation parameterized by 

Lbol = L44 (^v)° ' {1 + Z)< !0 44/j " 2 erg sec- 1 . (3) 

We take L44 = 2.9, a = 3.4 and ( = as a fiducial set of parameters on the 
basis of recent observational indications (David et al. 1993; Ebeling et al. 1996; 
Ponman et al. 1996; Mushotzky & Scharf 1997). Then we translate L bol (T gas ) 
into the band-limited luminosity L band [T gas , E\, E 2 ] by properly taking account 
of metal line emissions (Masai 1984) in addition to the thermal bremsstrahlung 
(we fix the abundance of intracluster gas as 0.3 times the solar value). 

Assuming that the intracluster gas is isothermal, its temperature T gas is 
related to the total mass M by 



fim p GM 



3r vir (M,z / )' 



where k B is the Boltzmann constant, G is the gravitational constant, m p is the 
proton mass, (i is the mean molecular weight (we adopt /u = 0.59), and 7 is 
a fudge factor of order unity which may be calibrated from hydrodynamical 
simulations or observations. The virial radius r v ; r (M, zj) of a cluster of mass M 
virialized at Zf is computed from A v i r , the ratio of the mean cluster density to the 
mean density of the universe at that epoch. We evaluate this quantity using the 
formulae for the spherical collapse model presented in Kitayama & Suto (1996b) 
and assuming for simplicity that Zf is equal to the epoch z at which the cluster 
is observed. Finally, we compute the mass function um(M, z)dM in equation 
(U) using the Press-Schechter theory (Press & Schechter 1974) assuming z = Zf 
as above. The effect of zj ^ z is discussed by Kitayama & Suto (1997) in this 
context, and the more general consideration of Zf 7^ z is given in Lacey & Cole 
(1993), Sasaki (1994), and Kitayama & Suto (1996a,b). 



3. Breaking the degeneracy 

The X-ray LogA^-LogS" for various Cold Dark Matter (CDM) models and the 
resulting constraints on S7o and erg are summarized in Fig|l|. Figure p] presents 
a clear example of the degeneracy among viable cosmological models in a sense 
that for a given value of Qq, one can find a value of as which accounts for the 
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Figure 2. The LogA-LogS relations of galaxy clusters for CDM 
models in (a) the soft X-ray (0.5-2.0 keV) band, (b) the hard X-ray 
(2-10 keV) band, and (c) the submm (0.85 mm) band. Lines repre- 
sent the models listed in table 1; L03 (solid), O045 (dotted), El (short 
dashed), L037 (long dashed), and LOla (dot-dashed). Also plotted in 
panel (a) are the la error bars from the RDCS (Rosati et al. 1995, 
1997), and the error box from the BCS (Ebeling et al. 1997). In panel 
(b), the arrows indicate the LogA-Log<S of all X-ray sources in the 
2-10 keV band from ASCA (Cagnoni et al. 1997), and the error bar 
(1 a) is the number of clusters observed by HE AO 1 (Piccinotti et al. 
1982). 

X-ray cluster LogA-LogS" . Several examples of such models are listed in Table 
1. 

Table 1. CDM models from the ROSAT X-ray Log A - Log 5. 



Model 


ft 


Ao 


h 


0-8 


a 


7 


L03 


0.3 


0.7 


0.7 


1.04 


3.4 


1.2 


O045 


0.45 





0.7 


0.83 


3.4 


1.2 


El 


1.0 





0.5 


0.56 


3.4 


1.2 


L037 


0.3 


0.7 


0.7 


0.90 


3.4 


1.5 


LOla 


0.1 


0.9 


0.7 


1.47 


2.7 


1.2 



The degeneracy among viable cosmological models can be broken by observ- 
ing wider (i.e., increase the statistics using wide-field surveys like SDSS, 2dF, 
PLANCK) and deeper (at higher redshifts) in different bands, (e.g., Eke, Cole, 
& Frenk 1996; Barbosa, Bartlett, Blanchard, & Oukbir 1996; Fan, Bahcall, & 
Cen 1997; Kitayama, Sasaki, & Suto 1998). Figure |2| plots the cluster LogA- 
Log5 predictions in soft-Xray, hard X-ray, and submm (due to the Sunyaev- 
Zel'dovich effect) for models which reproduce the observed LogA-LogS 1 in the 
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Figure 3. Redshift evolution of the number of galaxy clusters. Upper 
panels show the cumulative number iV(> 5, < z) against z in (a) the 
soft X-ray (0.5-2 keV) band with S = 10~ 13 erg cm -2 s _1 , (b) the 
hard X-ray (2-10 keV) band with S = 10 -13 erg cm~ 2 s _1 , and (c) 
the submm (0.85 mm) band with S u = 50 mJy Lower panels (d)-(f) 
are similar to (a)-(c) except for plotting the differential distribution 
dN(> S, < z)/dln z normalized by iV(> S). 



soft-Xray band (Table 1). Figure || illustrates the extent to which one can break 
the degeneracy between as and in CDM models {n = 1, h = 0.7) using the 
multi-band observational data. 

Similarly the redshift-distribution of cluster abundances can be a very pow- 
erful discriminator of the different cosmological models. Figure || exhibits the 
redshift evolution of the number of clusters in different bands. As expected, the 
evolutionary behavior strongly depends on the values of f2o and a$; the frac- 
tion of low redshift clusters becomes larger for greater Qq and smaller ag. It is 
indicated that one may be able to distinguish among these models merely by 
determining the redshifts of clusters up to z ~ 0.2. For this purpose, we adopt 
tentatively the X-ray brightest Abell-type clusters (XBACs, Ebeling et al. 1996), 
which is about 80 % complete and consists of 242 clusters with S > 5 x 10 -12 erg 
cm~ 2 s in the 0.1-2.4 keV band and z < 0.2. Keeping in mind the incomplete- 
ness of the sample and uncertainties especially in the estimated temperature 
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Figure 4. Contour maps on the ^lo-^s plane in (a) spatially flat 
(A = 1 - fio) CDM models, (b) open (A = 0) CDM models, (c) spa- 
tially flat CDM-like models with the fixed shape parameter (r = 0.25), 
and (d) open CDM-like models with T = 0.25. In all cases, h = 0.7, 
a = 3.4, and 7 = 1.2 are assumed. Shaded regions represent the la 
significance contours derived in KS97 from the soft X-ray (0.5-2 keV) 
LogA r -Log5 . Dotted and solid lines indicate the contours of the num- 
ber of clusters greater than 5 per steradian (10 2 , 10 3 , 10 4 from bottom 
to top) with 5 = 10 -13 erg cnr 2 s in the hard X-ray (2-10 keV) 
band and with S u = 50mJy in the submm (0.85 mm) band, respec- 
tively. Thick dashed lines represent the COBE 4 year result computed 
from the fitting formulae at 0.2 < U < 1 by Bunn & White (1997). 



data, however, we simply intend to perform a crude comparison with our pre- 
dictions in Fig. [|. 

The sky coverage of the XBACs is hard to quantify mainly due to the 
uncertain volume incompleteness of the underlying optical catalogue as noted by 
Ebeling et al. (1996). In Fig. ||, therefore, we simply plot the real numbers of the 
XBACs and normalize all our model predictions to match the total number of the 
XBACs at its flux limit 5(0.1-2.4 keV) = 5x 10 -12 erg cm -2 s _1 and redshift limit 
z = 0.2. In general, the model predictions are shown to be capable of reproducing 
well the shape and amplitude of the observed distributions N(> 5, > T, < z) 
even when T and z are varied. Taking into account the incompleteness of the 
observed data and large statistical fluctuations at low numbers, the agreements 
with models L03 and LO37 (both has = 0.3) are rather remarkable. Since 
the shapes and amplitudes of the predicted curves in these figures are primarily 
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determined by the value of Qq? this result provides a further indication for low 
fio universe. 




Figure 5. Tentative comparison with the XBACs (Ebeling et al. 
1996). Upper panel (a) shows the LogiV-LogS' in the 0.1-2.4 keV band 
with the redshift limit of z = 0.2 (models L03 and O045 almost overlap 
with L037 and El respectively). Middle panels exhibit iV(> S,>T,< z) 
versus z for 5(0.1-2.4 keV) = 5 x 10~ 12 erg cm~ 2 s _1 in the cases of 
(b) T = 2keV and (c) T = 5keV. Lower panels plot N{> S,>T,< z) 
versus T for 5(0.1-2.4 keV) = 5 x 10~ 12 erg cm~ 2 s~ 1 in the cases of 
(d) z = 0.2 and (e) z = 0.1. All the model predictions are normal- 
ized to reproduce the total number of the XBACs at its flux limit 
5(0.1-2.4 keV) = 5 x 10~ 12 erg cm" 2 s _1 and redshift limit z = 0.2. 
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4. Summary 



Let us summarize the conclusions of the present talk. 

(1) There exist several theoretical models which successfully reproduce the 
observed LogA^-LogS* relation of galaxy clusters over almost four orders 
of magnitude in X-ray flux. 

(2) The resulting as is given by the following empirical fit (95% confidence 
limit): 

a s = (0.54 ± 0.02 ± 0.1) x ^-o^-o.soQo+o.^ (g) 
for A = 1 - n CDM, and 

fjg = (0.54 ± 0.02 ± 0.1) x ^ o -0-^-o.9in 0+ o.68^ (6) 
for A = CDM. 

(3) Low-density CDM models (n = 1) with (n , A , h, a 8 ) = (0.3,0.7,0.7,1) 
and (0.45,0,0.7,0.8) simultaneously account for the cluster LogiV-LogS' , 
XTF, the COBE 4 year normalization. 

Maybe the most important point is that many cosmological models are more 
or less successful in reproducing the structure at redshift z ~ by construction. 
This is because the models have still several degrees of freedom or cosmological 
parameters which can be appropriately adjusted to the observations at z ~ (f2o, 
as, h, Xq, b(r,z)). We have shown that, given a complete flux limited cluster 
sample with redshift and/or temperature information, one can further constrain 
the cosmological models. In fact, our tentative comparison indicates that our 
predictions reproduce well the evolutionary features of the XBACs and that the 
results, although preliminary, seem to favor low density (0,q ~ 0.3) universes. 
As indicated by this preliminary result, surveys of objects at at high redshifts 
in several different bands (X-ray, radio and submm) are the most efficient and 
rewarding to break the degeneracy among the viable cosmological models. 
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